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Charge form factors corresponding to proton density distributions in exotic nuclei, such as 
6,8 He, 11 Li, 17 > 19 B and 14 Be are calculated and compared. The results can be used as tests of 
various theoretical models for the exotic nuclei structure in possible future experiments using a 
colliding electron-exotic nucleus storage ring. The result of such a comparison would show the 
effect of the neutron halo or skin on the proton distributions in exotic nuclei. 



1 Introduction 

It has been shown (see, e.g. the review [1]) that the study of exotic nuclei is likely to reveal new 
aspects of nuclear structure. Since the first experiments with radioactive nuclear beams (RNB) 
took place (e.g. [2]) it has been found from analyses of the total interaction cross sections that 
weakly-bound nuclei have increased sizes that deviate substantially from the R ~ 4~ x / 3 rule. 
It became clear (e.g. [3, 4]) that such new phenomena are due to the weak binding of the last 
few nucleons which form a diffuse nuclear cloud due to quantum-mechanical penetration (the 
so called nuclear halo). Another effect is that the nucleons can form a neutron skin [5] when 
the neutrons are on average less bound than the protons. Most of the neutron haloes have been 
studied for the lightest elements, such as He, Li, Be, B and other isotopes. 

Direct reactions are a well-known tool to provide information on nuclear structure. Most 
exotic nuclei are so shortlived that they cannot be used as targets. Instead, direct reactions with 
RNB can be done in inverse kinematics where the role of beam and target are interchanged. For 
example, the proton elastic-scattering angular distributions were measured for several incident 
energies, namely 25.1A MeV [6, 7], 41.64 MeV [8, 9, 10] and 714 MeV [11, 12] for 6 He, 26.254 
MeV (the preliminary data obtained in Dubna [13]), 324 MeV [11, 12], 664 MeV [11, 12] 
and 734 MeV [11, 12, 14, 15] for 8 He, as well as for energies less than 1004 MeV [11, 1] in 
the cases of 6 Li, n Li and others. In some of the theoretical analyses (e.g. [11, 12, 14, 16]) 
the eikonal approach using neutron and proton density distributions as well as parametrized 
nucleon-nucleon total cross sections have been used. The real part of the optical potential for 
calculations of 6 He+p, 6 He+ 4 He (£ Za6 =151 MeV) [17] and 6 He+p, 8 He+p (E inc < 1004 MeV) 



1 



[18] elastic scattering angular distributions was obtained microscopically using the realistic M3Y- 
Paris effective interaction [19, 20, 21] as well as the Tanihata neutron and proton densities of 
helium isotopes [22] in [17, 18] and the densities derived from the cluster-orbital shell model 
approximation (COSMA) [23, 11, 12, 24] in [18]. It was shown that elastic scattering is a 
good tool to distinguish between different density distributions [18]. The results of [18] were 
compared also with those from the alpha-core approach with the complex and fully non-local 
effective interaction [25] and the no-core model based on the large space shell-model calculations 
(LSSM)([26, 9, 27] and refs. therein). For larger energies(> 500 MeV/nucleon) it is accepted 
that the Glauber approximation is a suitable method to study charge and matter distributions 
from proton elastic scattering data [1, 28]. 

As mentioned in [1], however, an accurate determination of charge distributions of exotic 
nuclei can be obtained from electron-nucleus scattering in inverse kinematics using a colliding 
electron-exotic nucleus storage ring. 

The aim of the present paper is to calculate the charge form factors of exotic nuclei such 
as 6 He, 8 He, n Li, 17 B, 19 B and 14 Be on the basis of charge density distributions calculated by 
using various theoretical models. This work can serve as a challenge for future experimental 
works on measurements of the charge form factors and thus for accurate determination of the 
charge distributions in these nuclei. The latter can be a test of the different theoretical models 
used for obtaining charge distributions. Second, this work can throw light on the question about 
the effect of the neutron halo or skin on the proton distributions in the helium, litium, boron 
and beryllium isotopes. 

A brief representation of the theoretical scheme is given in Sec. II. The results and discussion 
are presented in Sec. III. 



2 The theoretical scheme 

The charge form factor F c /j(q) is the Fourier transform of the charge density distribution p c h( r ) : 

F ch ( q ) = [ dr Pch (r)e^ 1 (1) 
J J Pch(r)dr 

Let us firstly consider p c h( r ) obtained from a total nuclear wave function in which the positions 
r of the nucleons are defined with respect to the centre-of-mass of the nucleons in the nucleus 
(cm.). If the charge density is spherically symmetrical and normalized to unity, then: 

F ch (q) = Air J p c h(r)j (qr)r 2 dr, (2) 

where jo(qr) is the zero-order spherical Bessel function. Neglecting the charge neutron form 
factor, the charge density p c h{ r ) is represented by the density distribution of point-like protons 
in the nucleus p p0 mt(r) folded with the charge distribution of the proton /9 pr (r): 

Pch(r) = J ppointir - r')p pr (r')dr'. (3) 

Then the charge form factor (2) can be obtained in the form: 

F ch (q) = F point (q)F pr (q), (4) 



where F point (q) and F pr (q) are the form factors which correspond to the densities p P oint 

(r) and 

Ppr(r) (in formulae such as Eq. (2)), respectively. Secondly, in some calculations the wave 
functions which lead to the density distributions are often translationally non-invariant and the 
positions of the nucleons are determined with respect to the centre of the potential related to the 
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laboratory system (lab.) but not to the centre-of-mass of the nucleons in the nucleus. These wave 
functions admit a motion of the centre-of-mass which has to be switched off in the calculations. 
It is shown (e.g. in [29, 30, 31]) that the relationship between the form factor calculated by means 
of harmonic-oscillator wave functions (with all the nucleons in the Is state) K, mt,Zafe( c l) an d the 
form factor calculated by actual wave functions (with positions of the nucleons determined in 
the cm. system) -F po j ntiC m .(q) contains a correction factor: 

Fpoint,c.m.{c\) = / 4A F pointt i ab (q), (5) 
where b is the harmonic-oscillator parameter, or 

Fpoint, c. rn. (q) = 6^ ^ ^ Fpoint,lab(S\) i (6) 

R being the root-mean-square (rms) radius of the nucleus. For shell-model potentials different 
from the harmonic-oscillator one the Eqs. (5) and (6) are approximate. Thus in this case the 
expression for the form factor which can be compared with the experimental form factor can be 
written in the form [31]: 

FM = e^ 2 l &A F pomtMb {^)F P M). (7) 
In the Eqs. (1), (2), (4) and (7) all form factors are equal to unity at q=0. 



3 Results and discussion 

In this Section we calculate F po i nt (q) using the point charge density p po i n t(r) in Eq. (2) and 
-F c /i(q) using Eq. (4) for the exotic nuclei 6 ' 8 He, 11 Li, 14 Be and 17,19 B. For comparison we present 
also the charge form factors of the "non-exotic" 4 He and 6 Li nuclei. 

Various charge distributions of the proton p pr ( r ) can be used in order to obtain the form 
factor F pr (q). We use the Gaussian form normalized to unity: 

PpriT) = (7ra 2 )3 /2 exp(-r 2 /a 2 ) (8) 
which leads to the proton form factor 

F P M) = e~ q2aV4 . (9) 

For the value of the parameter a we choose a = 0.6532 fm which corresponds to the proton rms 
radius R p = 0.80 fm [32]. 



3.1 Density distributions of 4 ' 6,8 He 

We use firstly the following point nucleon densities p po int{f) for 4 ' 6,8 He deduced by Tanihata et 
al. [22] which determine F po i nt (q) in Eq. (4): 
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where X = Z, N. These densities were obtained by fitting the measured total reaction cross 
sections of 4 ' 6,8 He at 800^4 MeV incident on the C target using the optical limit of the Glauber 
model [33]. The parameter values a and b are determined from: 
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where a* = 1.53 fm for all 4 < 6 > 8 He isotopes, b* = 2.24 fm for 6 He and b* = 2.06 fm for 8 He [22]. 
Hence, a = 1.40 fm and b = 2.04 fm for 6 He, a = 1.43 fm and b = 1.93 fm for 8 He and a = 1.325 
fm for 4 He. Thus, the rms radii of the point-proton density distributions are 1.62 fm, 1.72 fm 
and 1.76 fm for 4 He, 6 He and 8 He, respectively. 

We use secondly the COSMA point nucleon densities [24] of 6 ' 8 He which have the same 
analytical form (10) and are based on the assumption of the 1^3/2 state for the relative motion 
of the alpha-core and each of the valence neutrons. The parameter values are [11, 12] a = 1.55 
fm, b = 2.24 fm for 6 He and a = 1.38 fm, b = 1.99 fm for 8 He. The rms radii of the point-proton 
density distributions are 1.89 fm and 1.69 fm for 6 He and 8 He, respectively. Here we would like to 
note the different trend in the behavior of the rms radii of the point-proton density distributions 
in the Tanihata density [22] and COSMA density [24, 11, 12]. While in the former the rms 
radius increases from 4 He to 6 He and 8 He, in the latter it decreases from 6 He to 8 He. Thus 
both densities contain opposite effects of the additional neutrons on the charge distributions in 
the helium isotopes. We note also that both Tanihata and COSMA densities have a Gaussian 
asymptotic behavior which is supposed not to be a realistic one for these nuclei at high q. Both 
densities had been used to fit mainly the total reaction cross sections and the rms radii. That is 
why it is necessary to consider also other more realistic theoretical predictions for the densities 
of the helium isotopes. We use also the proton densities of 4 ' 6 ' 8 He obtained within the LSSM 
calculations in a complete Ahuj shell-model space [26]. In them Woods-Saxon single-particle 
(s.p.) wave functions have been used for 6 He and 8 He and harmonic-oscillator ones for 4 He. We 
use also the experimental charge density for 4 He [34] with charge rms radius equal to 1.696(14) 
fm. 

3.2 Density distributions of 6 Li and n Li 

The COSMA point nucleon densities used in our calculations for 6 Li and n Li can be written in 
the following form [11]: 



x AT exp(-r 2 /a 2 ) 2 exp(-r 2 /6 2 ) [ „ 2 B (i ' 

P*oint{r) = N cX + N vX 2^ Ar* + - (r> - -6 2 



, X = Z,N. (12) 



For 6 Li (in the a+2N model): A=l, B=0, N cZ = N cN = 2, N uN = N vZ = 1, a=1.55 fm and 
6=2.07 fm [12]. For n Li (in the 9 Li+2A model): ,4=0.81, 5=0.19, N cZ = 3, N cN = 6, N vZ = 0, 
N uN = 2, a=1.89 fm and 6=3.68 fm [11]. The rms radii of protons are 2.44 fm and 2.31 fm 
for 6 Li and n Li, respectively. Similarly to the case of the helium isotopes, the COSMA density 
predicts a decrease of the proton rms radius from 6 Li to 11 Li when the number of neutrons 
increases. Secondly, we use the proton density of 6 Li obtained within a LSSM in a complete 
4huj shell-model space and of n Li in a complete 2hu> shell-model calculations [35]. For n Li 
Woods-Saxon s.p. wave functions have been used and harmonic-oscillator ones for 6 Li. Third, 
for 6 Li we use as well the density distribution of point-like protons [36, 37] which leads to the 
experimental charge distribution with rms radius equal to 2.57 fm [36]. 

3.3 Density distributions of 1719 B and 14 Be 

The point nucleon distributions in 17 ' 19 B and 14 Be were determined in [38] from analyses of the 
interaction cross sections of light radioactive nuclei close to the neutron drip line measured at 
around 800A MeV using the optical limit to Glauber type calculations by employing a harmonic- 
oscillator type densities including the contribution from the sci-shell: 
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where 
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(14) 



(15) 



and A is the width parameter. Its values are determined by a fit to the experimental interaction 
cross sections [38, 2]. In our work we use the values of A which correspond to the values of the 
effective rms radii of the point nucleon distributions obtained in [38]. They are 1.9366 fm, 2.0365 
fm and 2.0218 fm for 17 B, 19 B and 14 Be, respectively. The rms radii of the proton distributions 
are correspondingly 2.722 fm, 2.872 fm and 2.755 fm for 17 B, 19 B and 14 Be. Similarly to the 
case of the Tanihata densities of the helium isotopes the proton rms radius increases from 17 B 
to 19 B. 



3.4 Results of calculations 

In this subsection we present and compare the charge form factors F c h{q) for different nuclei 
calculated using Eq. (4) (which accounts for the proton size and in which the positions of the 
protons are determined from the centre-of-mass of the nucleons in the nucleus) and also make a 
comparison in one case with a form factor for point-like protons. 

In Fig. 1 are compared the charge form factors F c h(q) (4) of 6 He and 8 He calculated using 
the Tanihata [22] and the COSMA [11] densities. For comparison are given the charge form 
factors for 4 He calculated using the experimental charge distribution from [34] and the Tanihata 
density [22]. In Fig. 2 the charge form factors for 4,6 ' 8 He calculated using the LSSM densities 
[26] are compared with those for 4 He calculated with the experimental p c h{ r ) [34] and with the 
Tanihata density [22]. In Fig. 3 the charge form factors from (4) for 6 Li and n Li calculated 
using the COSMA densities [11, 12] are given together with those from calculations using the 
LSSM proton densities [35]. For comparison are given the charge form factors for 6 Li calculated 
using the experimental charge distribution from [36, 37] as well as the point-proton form factor 
for n Li calculated with the COSMA density distribution. The charge form factors of 17 B and 
19 B are compared in Fig. 4 with that of 14 Be using Eq. (4) and the point-proton densities from 
[38]. The charge rms radii and the diffuseness parameters of the charge densities are given in 
Table 1. The diffuseness parameter is defined as the distance over which the value of p c h{f) 
decreases from 90% to 10% of the value of the density in the centre of the nucleus p c h{ r = 0) 
divided by 4.4. 



3.5 Discussion and conclusions 

The main features of the obtained charge form factors can be summarized as follows: 

Although the number of the protons within a given group of isotopes is the same, the 
charge form factors are quite different from each other. This concerns 4 ' 6 ' 8 He nuclei (using the 
Tanihata, the COSMA and the LSSM densities), 6 Li and n Li (using the COSMA and the LSSM 
densities), as well as 17 B and 19 B isotopes (using the Suzuki densities). These differences reflect 
the differences in the phenomenologically obtained proton density distributions (which have 
been fitted to the empirical data for the total reaction cross sections) and in the corresponding 
radii. Thus the charge form factors can show the effects of the different number of neutrons 
and their distribution in the nucleus on the charge distribution in a given group of isotopes. 
Indeed, e.g. 4 He is a strongly bound a-cluster, while 6 He is a loosely bound a+n+n system. 
6 Li is predominantly an a+d formation while n Li is supposed to be a 9 Li+n+n system. We 
should emphasize, however, the different trend of the behavior of the charge form factors for 
6 He and 8 He when using the Tanihata, the COSMA and the LSSM densities. This follows, as 
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l < 1 < r 




Figure 1: Charge form factors from Eq. (4): a) for 4 He calculated with the experimental p c h( r ) 
[34] (dotted line) and with the Tanihata density [22] (dash- two-dotted line); b) for 6 He calculated 
using the Tanihata [22] (thick solid line) and the COSMA [11] (thick dashed line) densities, and 
c) for 8 He calculated using the Tanihata [22] (thin solid line) and the COSMA [11] (thin dashed 
line) densities. 

can be expected, the changes of the proton rms radii (i.e. the changes of the neutron effect on 
the charge distributions) within a given theoretical approach. While in the case of the Tanihata 
density the proton rms increases from 1.72 fm for 6 He to 1.76 fm for 8 He, in the COSMA and 
LSSM cases it decreases from 1.89 fm to 1.69 fm (in COSMA) and from 1.95 fm to 1.92 fm (in 
LSSM) for 6 He and 8 He. 

We should mention that the Tanihata, COSMA and Suzuki densities used in our work 
have Gaussian asymptotic behavior which is supposed not to be a realistic one for the nuclei 
considered at high q. This imposes the usage of other more realistic theoretical densities that 
may become available. We used the LSSM ones for 4,6 ' 8 He and 6 ' n Li nuclei. For 6 He, 8 He and 
n Li Woods-Saxon s.p. wave functions have been used in the LSSM calculations. 

We note also the important effect of the proton size accounted for in the calculations of the 
charge form factors. 

In summary, various existing models of charge densities for light exotic nuclei differ very 
much in their form factor predictions. They even predict opposite effects of the neutron excess 
on the charge distribution that could be tested by measuring charge form factors. Therefore, it 
would be most desirable to compare the form factors obtained here with the future experimental 
data on electron-exotic nuclei scattering. Our numerical results for the charge form factors are 
available upon demand. 
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Figure 2: Charge form factors for q > 1 fm^ 1 from Eq. (4): a) for 4 He calculated with the 
experimental p c h{ r ) [34] (dotted line), with the Tanihata density [22] (solid line) and using the 
LSSM density [26] (dashed line); b) for 6 He calculated using the LSSM density [26] (dash-dotted 
line), and c) for 8 He calculated using the LSSM density [26] (dash- two-dotted line). 
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Table 1: Charge radii (in fm) and diffuseness (in fm) of the charge densities. 



Nuclei 


Charge rms 


Diffuseness 


Densities 


4 He 


1.68 


0.30 


Exp. [34] 




1.81 


0.41 


Tanihata [22] 




2.03 


0.44 


LSSM [26] 


6 He 


1.89 


0.42 


Tanihata [22] 




2.06 


0.45 


COSMA [11] 




2.08 


0.49 


LSSM [26] 


8 He 


1.93 


0.43 


Tanihata [22] 




1.87 


0.41 


COSMA [11] 




2.04 


0.45 


LSSM [26] 


6 Li 


2.57 


0.54 


Exp. [36, 37] 




2.57 


0.52 


COSMA [11] 




2.63 


0.55 


LSSM [35] 


n Li 


2.45 


0.54 


COSMA [11] 




2.51 


0.57 


LSSM [35] 


14 Be 


2.87 


0.64 


Suzuki [38] 


17 B 


2.84 


0.60 


Suzuki [38] 


19 B 


2.98 


0.63 


Suzuki [38] 
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